arXiv:1503.06615v2 [cond-mat.soft] 15 Jun 2015 


The effective colloid interaction in the Asakura-Oosawa model. Assessment 
of non-pairwise terms from the virial expansion 

Andres Santos,^'Mariano Lopez de Haro,^’'^^ Giacomo Fiumara,^'and Franz Saija^’'^^ 

Departamento de Fisica and Instituto de Computacion Cientifica Avanzada (ICCAEx), Universidad de Extremadura, 
E-06071 Badajoz, Spain 

Instituto de Energias Renovables, Universidad Nacional Autonoma de Mexico (U.N.A.M.), Temixco, 

Morelos 62580, Mexico 

Department of Mathematics and Computer Science, University of Messina, Viale F. Stagno DAlcontres 31, 
1-98166 Messina, Italy 

CNR-IPCF, Viale F. Stagno d’Alcontres, 37-98158 Messina, Italy 
(Dated: 8 August 2016) 

The relevance of neglecting three- and four-body interactions in the coarse-grained version of the Asakura- 
Oosawa model is examined. A mapping between the first few virial coefficients of the binary nonadditive hard- 
sphere mixture representative of this model and those arising from the coarse-grained (pairwise) depletion 
potential approximation allows for a quantitative evaluation of the effect of such interactions. This turns out 
to be especially important for large size ratios and large reservoir polymer packing fractions. 


I. INTRODUCTION 

The description of the thermodynamic properties of 
complex fluids is in general not an easy task. In it, one 
faces the presence of many degrees of freedom and maybe 
also of different length and time scales. An approach that 
is usually followed is to recur to coarse-graining. In this 
approach, what one attempts is to integrate out the irrel¬ 
evant degrees of freedom and hence to end up with a sim¬ 
pler (equivalent) system with an effective interaction that 
hopefully captures exactly the essential features of the 
real interaction. Achieving an exact coarse-graining is, 
however, also difficult. This is due to the fact that, even 
if the underlying original molecular interactions are pair¬ 
wise, the resulting effective potential turns out to be in 
general a many-body one. For this reason, in the coarse- 
graining process it is usual to replace the full many-body 
potential by a simpler effective one in which only pair 
interactions are involved. The question then arises as to 
whether the thermodynamic properties derived with the 
effective potential provide a reliable account of the same 
properties for the original fluid. 

In the early 20th century, Kamerlingh Onnes^ intro¬ 
duced the (then empirical) virial series to provide a math¬ 
ematical representation of experimental pressure-density- 
temperature data of gases and liquids. In a broader con¬ 
text, one refers to a virial expansion of a given property 
when such a property is expressed as a power series in 
density. Thus, the virial expansion represents in princi¬ 
ple a systematic way for calculating the properties of bulk 
matter, provided of course that the coefficients in the ex- 
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pansion (the so-called virial coefficients) are known accu¬ 
rately either through direct measurement or from theo¬ 
retical developments. At least formally, one of the great 
achievements of statistical mechanics, and a major break¬ 
through in the theoretical approach to calculating virial 
coefficients involved in the equation of state of fluids, 
occurred when Mayer^’^ was able to obtain general ex¬ 
pressions for the corresponding virial coefficients in terms 
of sums of cluster integrals over the interaction among 
groups of fluid particles. In particular, the second virial 
coefficient depends only on pair interactions, while the 
third virial coefficient depends on two- and three-body 
interactions, and so on. If the intermolecular potential 
is simple enough, some of such virial coefficients may be 
calculated analytically. For instance, in the case of the 
hard-sphere (HS) fluid, the first four virial coefficients 
are known analytically. The same applies to the cele¬ 
brated Asakura-Oosawa (AO) model,which describes 
colloidal HSs in a solvent of ideal polymers that inter¬ 
penetrate each other but interact with the colloids via 
a HS repulsion. In this case, considering that the sys¬ 
tem may be taken as a nonadditive hard-sphere (NAHS) 
binary mixture with high size asymmetry, analytical re¬ 
sults for the first four virial coefficients have been very 
recently reported.^ In general, however, numerical eval¬ 
uation is required and even for simple potentials such 
as the ones just mentioned there are various technical 
difficulties involved in computing the higher order virial 
coefficients. A noteworthy aspect of the usefulness of 
virial coefficients is that the comparison of these coeffi¬ 
cients for the original system and the ones stemming out 
of the coarse-grained interaction potential will indicate 
to what extent and under which conditions the neglect 
of the many-body terms in the latter has an impact on 
whether the thermodynamic properties of both systems 
agree. 

Although largely ignored for about 20 years, interest in 
the AO model grew in the 1970s and 1980s and it started 
to get significant attention in the 1990s, which contin- 
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ues up to this day 7“^® It is well known that its coarse¬ 
grained description involves an effective (depletion) pair¬ 
wise interaction between the colloids that, among other 
things, leads to fluid-fluid demixing. Attempts to exam¬ 
ine whether the thermodynamics obtained with the de¬ 
pletion potential agrees with the one of the full mixture 
have also been reported. In particular, it was found that, 
for a polymer/colloid size ratio q smaller than the thresh¬ 
old value go = 1 ~ 0.1547, the AO pair depletion 

potential turns out to be the only one contributing to 
the exact effective interaction among the solutes.^® Re¬ 
cently, Ashton and Wilding^®’^’’ focused, via simulation, 
on the dominant many-body effect neglected in the pair 
potential description if g > go, namely the one associated 
with the interaction between three colloidal particles. To 
this end, they examined the difference between the third 
virial coefficient of the full system and that of the effec¬ 
tive system. 

Here we will follow a similar route but also profit from 
the availability of the first five virial coefficients of the 
full NAHS mixture corresponding to the AO modeF to 
derive, by the exact mapping that may be performed be¬ 
tween both sets of coefficients, partial contributions to 
the second, third, and fourth virial coefficients of the ef¬ 
fective one-component colloidal fluid. The aim is to quan¬ 
tify the deviations from the exact results that one gets 
for both the third and fourth virial coefficients when com¬ 
puting them with the coarse-grained depletion potential. 
As we will see, the influence of non-pairwise interactions 
on the third and fourth virial coefficients is rather small 
for g < 0.4 but becomes increasingly important for larger 
values of the size ratio. 

The paper is organized as follows. In Sec. II, and in 
order to make the paper self-contained, we recall the re¬ 
sults for the first five virial coefficients of the original AO 
binary mixture. This is followed in Sec. Ill by the intro¬ 
duction of the osmotic pressure of the colloidal system, 
which allows us to make the mapping between the virial 
coefficients of the mixture and those coming out of the 
effective one-component colloidal system. The analytical 
results that follow from the (coarse-grained) pair deple¬ 
tion potential approximation are presented in Sec. IV. 
Section V provides a comparison between the exact and 
the approximate results. The paper is closed in Sec. VI 
with further discussion and some concluding remarks. 


II. VIRIAL COEFFICIENTS OF THE AO MODEL 

Consider a binary fluid mixture of V = + Np 

spheres (colloids-|-polymers) in a volume V. The colloid 
and polymer mole fractions are Xc = Nc/N and Xp = 
Np/N = 1 — Xc, respectively. Analogously, the partial 
and total number densities are pc = Nc/V, Pp = Np/V, 
and p = Pc + Pp = N/V. The interactions are assumed to 
be of HS type. The distance of closest approach between 
spheres of species a and 7 , denoted by <7aj, is such that 
CTcc = Cc, Cpp = 0 , and (Jcp = 50 -c(l -|- g), with the size 


ratio g acting as the (positive) nonadditivity parameter. 
The colloid packing fraction is rjc = f For simplic¬ 
ity, from now on we choose ctc = 1 as the unit of length. 
This NAHS mixture defines the well known AO model. 
The usual virial expansion of the mixture reads 


I3a{pc, Pp) =Pc In (PcA®) -h pp In (ppA®) 

BnjXc, g) n 
^ n — 1 ’ 

n—2 


p 


( 2 . 1 ) 


l3p{Pc,Pp) = P+^Bn{Xc,q)p'^, (2.2) 

n—2 


where a is the free energy per unit volume, p is the pres¬ 
sure, j3 = l/ksT (ks being the Boltzmann constant and 
T being the absolute temperature), Aq, is the thermal 
de Broglie wavelength of species a, and the notation 
Bn{xc,q) makes it explicit that the virial coefficients de¬ 
pend only on the mole fraction Xc of the colloids and on 
the size ratio g. 

The second, third, fourth, and fifth virial coefficients 
of the AO model are given by^ 

B 2 {xc, q) = xlBii + 2XcXpBi2{q), (2.3) 


B3,{Xc, q) = xlCiii + ^xlxpCii 2 {q), (2.4) 


Bi{xc,q) = x'^^Diiii + 4 a;®XpHi 112(g) -h &xlxlDii22{q), 

(2.5) 

B^{xc,q) =xlEiiiii + bXcXpEiiii2{q) 

+ 10a;®a;pAiii22(g)+10a;2a;pAii222(g), (2.6) 

where all the composition-independent coefficients, ex¬ 
cept E 11112 and i?iii 22 j are exactly known as functions 
of g. 


Bii = 



Bi2 — 


TT jl + q)^ 
6 2 


(2.7) 


Cn. = (^)bo, C„. = (I)*' + ®‘' + ‘''’' + *‘' 

(i) 


( 2 . 8 ) 


Am = 


'ttnS /2707 219 V 2 4131 cos"® i' 
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707r 
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TABLE I. Numerical values of the partial coefficients -Eirir2 
and -E11122 for some values of the size ratio q. The error on 
the last significant figure is enclosed in parentheses. 


q 

511112 

5ii122 

0.05 

0.0267(6) 

« -5 X 10”“ 

0.10 

0.0437(4) 

-4.3(9) X 10”® 

0.15 

0.0666(8) 

-3.4(8) X 10”® 

0.20 

0.0955(9) 

-1.6(2) X 10® 

0.30 

0.177(1) 

-0.00147(4) 

0.40 

0.296(4) 

-0.0075(2) 

0.50 

0.457(2) 

-0.0276(5) 

0.56 

0.575(5) 

-0.0546(5) 

0.60 

0.666(3) 

-0.0831(9) 

0.70 

0.931(4) 

-0.216(2) 

0.80 

1.257(7) 

-0.505(4) 

0.90 

1.652(6) 

-1.085(4) 

1.00 

5iiiii 

-2.18(3) 


cal ensemble (iVp, Nc, V, T). On the other hand, in order 
to analyze the effective one-component colloidal fluid, it 
turns out to be convenient to consider the semi-grand- 
canonical ensemble (/ip, N,., V, T), where /ip is the chem¬ 
ical potential of the polymer component. In that ensem¬ 
ble, the pressure of the mixture can be written as^^ 

Pvi.Pci Zp) = Zp /3n(pcj’Zp), (3-1) 

where 

g/S/ip 

Zp = (3.2) 

is the polymer fugacity and n(pc, Zp) is the osmotic pres¬ 
sure that takes into account the (formally) exact effective 
colloid-colloid interactions mediated by the polymers. Its 
virial expansion is 




27g® 162q^ 81g® 


^1112 ~ 


/7^^ 


V6/ 

2807rL 


35 
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~W 

( 2 . 12 ) 


-486Pi(Q^-h9) 


^ (log® - 5ig^ -b 2iog2 -b 6976) 

9+^0 rczrS oozr6 


-P 2 (5g® - 28Q® 


-bl 29 g'‘- 124 g 2 -b 11378) 


Pi 


= (0 65 , &5 - 28.224512, 


(2.13) 

(2.14) 


P 11222 =- QY ^ (3240 + 7695(7 + 6780(72 + 2706(7® 

+492(7^ + 419®). (2.15) 

In Eq. (2.13), Q = + 6(7 — 1, Pi = tan“® g, and 

P 2 = tan”® [g /(9 + 1 )]. 

To our knowledge there are no analytical results for the 
composition-independent coefficients P 11112 and P 11122 
for general values of q. Therefore, we have computed 
them by a standard Monte Carlo (MC) numerical inte¬ 
gration procedure for a number of values of q in the range 
0.05 < 9 < 1. The results are displayed in Table I, which 
is more extensive than the equivalent table of Ref. 7. 

There exist approximate analytical theories, like the 
free volume (FV) theory,®® that account in closed form for 
the equation of state of the full AO mixture. For further 
use, Appendix A provides the approximate expressions of 
the first few virial coefficients arising from the FV theory. 


III. OSMOTIC PRESSURE AND VIRIAL COEFFICIENTS 
OF THE EFFECTIVE COLLOIDAL SYSTEM 


Equations (2.1) and (2.2) are expressed in the canoni- 


/3n(pc, Zp) = Pc + XI : (3-3) 

n^2 

where B‘^{zp,q) are the virial coefficients of the effec¬ 
tive one-component colloidal fluid. The fugacity Zp of 
the polymer component can be equivalently represented 
by the reservoir polymer packing fraction rjp^^ = Zp^q^. 
Thus, henceforth we make the change B‘)^{zp, 9 ) —>■ 
B^{r]p^^,q). The effective virial coefficients B‘^{iqp^^,q) 
can be further expressed as a series in powers of 

00 

Bf{Vp,r,q) = X (3-4) 

j=0 

Our aim in this section is to provide the exact rela¬ 
tions between the effective one-component virial coeffi¬ 
cients Bn'^ (with n+j < 5) and the binary-mixture virial 
coefficients of Eqs. (2.3)-(2.6). The details are given in 
Appendix B with the results 

b!°'> =Bii, (3.5a) 

b’^^ =Diiii, Bi°^=Eiiiii, (3.5b) 

^ 2 '^ = ^ (^^®®2 - 2 ^ 12 ) , (3.6) 

(3.7) 

(^g £*1112 - 65120112 + 3 .^ 12 ^ j (3.8) 
B^^ = (55iii22 — 165i25ii22) , (3.9) 
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^ 4 ^^ — ^Bi2Diii2 - ^^112 

+ 9Bl^Cii2-2Bt2), (3.10) 

Equations (3.5)-(3.10) provide the sought relationships 
between the effective and binary-mixture virial coeffi¬ 
cients that account for all the three- and four-body in¬ 
teractions up to the fourth virial coefficients of the ef¬ 
fective system. Making use of Eqs. (2.7)-(2.15) and the 
numerical values of Table I, one can then know the ex- 
act g-dependence of the coefficients Bit . Furthermore, 
in the case of the approximate FV theory, the results are 
explicitly given by Eqs. (A6)-(A8). 


with fnsir) = -0(1 - r) and = 0(r - 1)0(1 -k 

q - r)[uj{r)]^ /j\. 

On the other hand, if go < 9 < 1; m-body terms with 
3 < m < 11 gradually contribute to $eff, the upper limit 
(m = 11) being due to the fact that a polymer can overlap 
simultaneously with 12 nonoverlapping colloids only if 
g > 1.^® Therefore, Eq. (4.2) becomes an approximation 
(henceforth referred to as the coarse-grained depletion 
approximation) if g > gg. In the remainder of this section 

( 1 ) 

we explicitly evaluate the effective virial coefficients Bit 
with n j < 5 for any 0 < g < 1 in this coarse-grained 
approximation. 


IV. VIRIAL COEFFICIENTS ARISING FROM THE USE 
OF THE EFFECTIVE PAIR AO POTENTIAL 


By integrating out the polymer degrees of freedom, it 
is possible to derive the formally exact effective many- 
body interaction potential of the colloids, <l>eff(r'^‘=) = 
<l’eff(ri, r 2 ,..., r^v,.), in the AO model. The result is^® 


r r 

=-Zp Mr 1-0 

Nc-1 Nc 

i=l j=i-il 


1 + 9 



(4.1) 


A. Second virial coefFicient 

As a consequence of (4.5), the second virial coefficient 


IS 


^OO 

Bf^ = -2TT drr^/dep(r) 

Jo 

O'TT 77 I 

= Y - 271 - ^ ^ drr^ [uj{r)Y 


From here one can easily obtain 


B. 


( 1 ) 


SHS 


= -^(1 + ^ + ^ + ^) 
^ 4 ^ 2 ^ 12 ' ’ 


(4.6) 


(4.7) 


where 0(a;) is the Heaviside step function and is 

the original colloid-colloid HS pair potential of diameter 
CTc = 1. If g < go, a polymer particle cannot overlap 
with more than two nonoverlapping colloids, so that 
is exactly given by^® 

,54>eff(r^'=) =-ZpV [1 - 77c(l -k g)^] 

Afc-l Me 

+ YY'depij'ij), (4.2) 

i —1 j—i +1 

where (/>dep(r) is the effective AO pair depletion potential. 
It is given by 


/3</'dep(7’) 


/ 

CX), 

< -7?p,rW(r), 


r < 1, 

1 < r < 1 -k g, 
r > 1 -k g, 


where 

w(r) = ^(1 -k g - r)2(2 -k 2g -k r). 


(4.3) 


(4.4) 


The corresponding Mayer function is 
/dep(r) =/Hs(r) + - l] 0(r - 1)0(1 + g - r) 

00 

pd) i^\^i ( 4 _ 5 ) 


=/hs( 7') + X! /hep(’^)^: 
7=1 


'p,r ’ 


Bf^ _ 27 

“ ~40g 


16g 68g2 17g® 17q^\ 

~9' ^ ^ ’ 

(4.8) 


277r 

“~ 112 g 2 V 



19g , 113g2 ^ 451g® , 41g4 
~8~ 54 540 

(4.9) 


where B^^ = Equations (4.7)-(4.9) agree with Eqs. 
(3.6) and (3.7). Of course, this is an expected result since 
the pair approximation is exact at the level of the second 
virial coefficient. 


B. Third virial coefFicient 


We now turn to the third virial coefficient 


3 q 


dk 


/dep(fc) 


where 


(4.10) 


/dep(fc) = J fdepir) 

Att 


dr r sin(fcr)/dep(7’) 


(4.11) 
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is the Fourier transform of fdepif)- From Eq. (4.5) we 
have 

OO 

Uep{k) = /Hs(fc) + ^ (4.12) 

i=i 


C. Fourth virial coefFicient 

Finally, we consider the effective fourth virial coeffi¬ 
cient in the coarse-grained pair depletion approximation. 
It is given by 


where 


~ 47r 

fmik) =-^{k cask-sink), (4.13) 


^ 1 47r /■1+® 

fdJpi^) drrsin(fcr) [w(r)]^ 


(4.14) 


According to Eqs. (4.10) and (4.12), the coefficients 
and are given by 


) = -(27r)- 


dk 


fmik) /^(fc), (4.15) 


Bf = - (2^)-3 


Hs(fc)S(fc) 


-b 


Sw]'}- 


(4.16) 


^ 4 ® = -^ (^3^ + 6 ^+ , (4.20) 

where each thick bond represents a Mayer function (4.5). 
Expanding in powers of one gets 

(4.21) 

Now, a thin solid line between two circles represents the 
HS Mayer function /hs(*'); while a dotted line represents 
a term f^Jp{r). Interestingly, Eq. (4.21) can be written 
as 


Bf ^ = -Gtt 


r-i-t? 


dr r^a;(r)j /2 (^)) 


(4.22) 


where 


The explicit expressions of f\]}pik) and f^]}pik) can be 
obtained from application of Eqs. (4.4) and (4.14) but, 
for conciseness, they will be omitted here. Insertion of 
those expressions into Eqs. (4.15) and (4.16) yields (for 
0 < 9 < 1 ) 

of. , ^<1 4g2 14g3 , 2q^ ^ q\ q^ \ 

BfS + 5 25 75 ^ 175 ^ 35 ^ 315^ ’ 

(4.17) 

bP 27 / 107q 1529g2 9253^3 3889^^ 

BfS “ “ ^ V “ “90 315 2520 5670 

4663g5 104996 10499^ \ 

18900 9450 113400/■ ’ 


is the HS cavity function to second order in density, which 
is exactly known. 

After some lengthy algebra it is possible to find a fully 
analytical expression for B^^ (see Appendix C). Again, 

since that expression of B^'^ is exact for 9 < 90 , use of 
Eq. (3.10) allows us to derive an exact analytical form of 
B 11112 for 9 < 90 [see Eq. (C9)]. The values correspond¬ 
ing to 9 = 0.05, 0.10, and 0.15 are 0.026 588 4,0.043 507 4, 
and 0.066 359 4, respectively. Comparison with the sec¬ 
ond column of Table I shows again an excellent agreement 
of the MC results with the exact values. 


As a byproduct, since Eq. (4.18) must be exact for 
9 < 9 o, Eq. (3.9) allows one to obtain the exact expression 
of B 11122 for q < qo, 

r. 7^496 f 799 28I92 352996 25199 ^ 

-^11122 — —- I H”-“1“-H”-H“- 

800 V 18 63 504 378 

3458396 375996 37699" \ 

11340 5670 68 040/ ■ ^ 

The exact values of B 11122 corresponding to 9 = 0.05, 
0.10, and 0.15 are -1.229 60 x 10"", -4.250 92 x 10 - 6 , 
and —3.500 57 x lO”®, respectively. By comparison with 
the third column of Table I, we observe that the MC 
results agree with the exact values within the associated 
uncertainties. 


V. COMPARISON BETWEEN THE EXACT AND THE 
APPROXIMATE COEFFICIENTS , AND Bf^'> 

As said before, coarse-grained pair-potential approxi¬ 
mation (4.2) is only correct if 9 < 90 ~ 0.1547. Beyond 
that value, the virial coefficients B®® (with n > 3) ob¬ 
tained from the AO pair potential (4.3) differ from the 
exact ones. In particular, the exact coefficient B|® is in¬ 
fluenced by three-body interactions, ^6 ,27 ^qg gxact 

coefficient B|® is influenced by both three- and four-body 
interactions. 

Here we restrict ourselves to Bg^^ and B^^^ In 
the coarse-grained approximation, they are given by Eqs. 
(4.17), (4.18), and (C4)-(C8), respectively (if 9 < 1). 
The exact expressions are given by Eqs. (3.8), (3.9), 
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FIG. 1. Comparison between the exact and the coarse-grained 
coefficients (a) (b) and (c) B^^K 


and (3.10), respectively, in terms of the composition- 
independent virial coefficients of the binary mixture. 
While in Eq. (3.8) all the coefficients are known ana¬ 
lytically, in Eqs. (3.9) and (3.10) one needs to resort (for 
q > go) to numerical MC evaluations listed in Table I. In 
the FV theory for the full AO binary mixture, and 
are given by Eqs. (A6) and (A7), respectively, while 

bP = 0 . 

The comparison between the exact and approximate 
coefficients is carried out in Fig. 1. We see that the influ¬ 
ence of three-body interactions on B^^^ and B^^ is prac¬ 
tically negligible in the range go < q ^ 0.4 but becomes 
quite important, especially in the case of B^^\ if 9 ^ 0.6. 



To aid visibility, the curves for q = 0.8 have been shifted 
vertically by 1.0. 


A similar conclusion can be drawn from B^^: the role 
played by three- and four-body interactions is irrelevant 
if <7 < 0.4 but becomes essential as q increases. We ob¬ 
serve that the non-pairwise contributions to the true ef¬ 
fective many-body colloid potential tend to increase the 
values of and i?|® with respect to the coarse-grained 
estimates, thus partially compensating for the attractive 
character of the pair depletion potential. For instance, 
while the coarse-grained approximation predicts a mono¬ 
tonic decrease of B^^'^ with increasing g, the exact coeffi¬ 
cient presents a non-monotonic behavior with a minimum 
at <7 ~ 0.54. Also, B^^'^ is negative definite in the coarse¬ 
grained approximation, while it actually changes from 
negative to positive at g ~ 0.74. As for the FV theory, 
it qualitatively agrees with the main trends of the exact 
coefficients B^^ and B ^^'^, especially as q increases. 

Assuming sufficiently small values of expansion 
(3.4) can be truncated to obtain the approximate forms 
BfiVp.r) ~ and Bf{iqp^,) ^ 


Bf^ + B^^^rjp^r- The resulting curves for q = 0.154 < go 
and g = 0.8 in the range 0 < ryp^r < 0.1 are plotted in 
Fig. 2. Figure 2(a) is qualitatively analogous to Fig. 3 of 
Ref. 27. We observe that the impact of three- and four- 
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FIG. 3. Plot of the osmotic compressibility factor piijpc ver¬ 
sus the colloid packing fraction r)c for = 0.1 and q = 0.4, 
0.6, 0.8, and 1.0. To aid visibility, the curves for q = 0.4, 0.6, 
and 0.8 have been shifted vertically by 0.05, 0.10, and 0.15, 
respectively. 


body interactions on and i?|® becomes relevant for 
Vp.i ^ 0-02 if g = 0.8. 

Notwithstanding the results displayed in Figs. 1 and 2, 
and since the coarse-grained description gives the exact 
effective second virial coefficient the osmotic pres¬ 
sure of the colloids in the presence of nonadsorbing poly¬ 
mers is expected to be well described by the coarse grain¬ 
ing method if the behavior is dominated by the second 
virial coefficient. To clarify this point, we consider the 
expansions (3.3) and (3.4) truncated for n-|-j > 6, which 
implies small values of both packing fractions rjc and rjp^v 
Figure 3 shows the corresponding compressibility factor 
PH/pc within the range 0 < rye < 0.1 for rjp ,. = 0.1 and 
several values of the size ratio q. Under those conditions, 
no difference between the exact and coarse-grained re¬ 
sults are visible for q = 0.4 and very small deviations 
can be observed for q = 0.6 near r]c = 0.1. Only for 
high size ratios {q = 0.8 and 1.0) it is apparent that 
the coarse-grained approximation underestimates the os¬ 
motic pressure; an effect that is expected to become more 
and more important as the packing fractions rjc and rjp^-c 
increase beyond the range of applicability of the trunca¬ 
tion for n -b j > 6. It is interesting to note that the FV 
theory is rather close to the exact results for q = 0.8 and 
1 .0. 


VI. CONCLUDING REMARKS 

Using the available results^ for the virial coefficients 
of the AO binary-mixture model, we have assessed, for 
size ratios 0 < g < 1, the effect of neglecting three- and 
four-body interactions on the values of the effective one- 
component virial coefficients and that follow 
from the depletion pair potential derived in the coarse¬ 


grained version of such a model. While it was already 
well known that the coarse-grained version is exact for 
q < qo = 2/-\/3 — 1 ~ 0.1547, the mapping between the 
virial coefficients of the true mixture and the effective 
ones that we have presented here, together with the corre¬ 
sponding analytical results, have allowed us to explicitly 
quantify the differences for the partial contributions , 

B^\ and B^^'^ for q > qo- As an extra bonus of this map¬ 
ping, exact analytical expressions for the binary-mixture 
coefficients E 11122 and E 11112 were derived for any size 
ratio q < qo- The same was in turn useful to check the 
accuracy of our numerical results for those coefficients, 
which were proven to be very reliable. 

The results indicate that the coarse-grained pair de¬ 
pletion approximation is very accurate for go < <7 ^ 0-4 
but one must certainly take into account the influence of 
three-body interactions on and i3|® if g > 0.6, their 
role becoming essential as g increases. While it is not pos¬ 
sible at this stage to disentangle the roles of three- and 
four-body interactions on i3|®, it is reasonable to expect 
that four-body terms could be important at least for val¬ 
ues of g close to unity. All these facts should be especially 
noteworthy when dealing with dense systems. Also, for 
small values of such an influence has been shown 
here to be relevant. In fact, as already pointed out in the 
case of three-body interactions by the numerical stud¬ 
ies of Ashton and Wilding,the deviation between 
the exact and coarse-grained values of significantly 
increases as 77p_r becomes larger. Therefore, care must 
be exercised when drawing conclusions from the coarse¬ 
grained version of the AO model if either g or r^p ,., or 
both, are large. 
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Appendix A: The FV theory 


In the FV theory^^’^"*’ the free energy of the system is 
expressed as a sum of a term corresponding to a pure 
colloidal suspension in the volume V and a term cor¬ 
responding to a pure polymer solution in the volume 
a{r]c)V, where the free volume fraction a(rjc) is motivated 
by scaled particle theory. The corresponding equation of 
state is^° 


P 


=XcZcs{Vc) + 
-|-3g(3 + 2g) 


Xpqi]c 


l-rjc (1 - rjcY L 
..2 


3 -|- 3g -f g 


Pc 


1 - Ifc 


9g^ 


Pc 


(1 - PcY 


2 

(Al) 












where Zcs(??) = (l+? 7 +? 7 ^—??^)/(l—is the Carnahan- 
Starling compressibility factor of a one-component HS 
fluid. Equation (Al) is consistent with the exact second 
and third virial coefficients [see Eqs. (2.3), (2.4), (2.7), 
and (2.8)]. On the other hand, the FV fourth and fifth 
virial coefficients are approximate only. They are given 
by Eqs. (2.5) and (2.6) with 


= (0'l8, £;iiiii = (0%8. 

(A2) 

/ TT \ 3 1-1-9g-1-36g^-1-30g® 

(A3) 

^ ,, /n\^ l + 12q + 66q^+ 76q^ 

Enn2{q)=[^) 5 

(A4) 

£> 1122 ( 9 ) = £' 11122 ( 9 ) = £ 11222 ( 9 ) = 0- 

Insertion into Eqs. (3.7)-(3.10) yields 

(A5) 

4'^ ( 9 ) =- (0 ' (24 + 33g - 3q^ - 20(7® - 

129 "^ 


(A6) 


where we have used Eqs. (2.3)-(2.6). Next, from Eq. 
(2.1) and the thermodynamic relation = {da/dpp)p^, 
we obtain 


3 4 

Ppp = In (ppAp) -|- 2 Bi 2 Pc + —Cii 2 Pc + gT^iii2Pc 

5 

+4:Dii22PcPp + ^Aiiii2P^ -I- 5Eiii22PcPp 
+ -^Eii222pipl 

Consequently, fugacity (3.2) can be written as 


~ —1 + 2Bi2Pc + ( 25^2 + aC'ii2 ) Pc ffi ( 


-|-3i?i2Cii2 -I- -7?1112j Pc + 4I?1122PcPp 
+ ^2^12 + 3^12^*112 + g.Bi2lllll2 + g^ll2 
+ ^Aiiii 2 ^ Pc + (8ill27?1122 + 5i5iii22) PcPp 


+ -^Aii222PcPp + 0{p^). 


(B3) 


=- (54 + ^ - 72g2 _ I47g3 - 

lUq^ 43g6 33g^ 3g® q^\ 

^ 8 ~Y~ 

^ 5(3) ^ ^(2) ^ 

Thus, only the coefficients B^\ and B^^ 

actly given by the FV theory. 


- 63g'‘ 

, (A7) 

(A8) 
are ex- 


Appendix B: Derivation of Eqs. (3.5)-(3.10) 

We start by rewriting Eq. (2.2) to fifth order as 

/3p(Pc) Pp) =Pc + Pp + Biip^ + 2Bi2PcPp + Ciiip® 

+ 3Cii2PcPp + T^llllPc + 4Zliii2PcPp 
+ 6Zlii22PcPp + Aiiiiip® -|- 5Eiiii2PcPp 
-I- 10Eiii22PcPp + 10Aii222PcPp + C>(p®), 

(Bl) 


This can be inverted to express pp as a series expansion 
in powers of Pc and Zp, 


^ =1 - 2Bi2Pc + ( 25^2 - -C 112 ] pi-( -Bt2 


— 3 B 12 C 112 + —Diii2j Pc — 4Z1ii22Pc^p 

+ ^ 3-^12 ~ ^Bf 2 Cll 2 + -B 12 D 1112 + gC '^12 

~^Aiiii 2 j Pc + (16i?1271ll22 — 5A11122) Pc^^p 


15. 


--^Eii222Pc2:p + ' 


(B4) 


where the ellipsis denotes terms of order PcZp with n+j > 

6 . 

Inserting Eqs. (Bl) and (B4) into Eq. (3.1) one can 
easily identify the coefficients shown in Eqs. (3.5)-(3.10). 
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Appendix C: Expressions for in the coarse-grained approximation 

In the range of interest l<r<l-|-g<2, the expressions for the contributions of y 2 ^{r) in Eq. (4.23) are^®’^° 


Vir) = 




tlj{r) = 


TT^ (r — 2)^ 
^ 35r 


(r® -b 4/ - Slr^ - lOr^ -b 479r - 8l) , 


(Cl) 

(C2) 




M 


=7r0 ( V 3 — r 




23 36 


_ J .2 _ ( -12L ] QOS ^ — — 

15 35r; y'3(4 - r^) 


3r® 2r 9 \ + r — 3 


”^'^560 15 2 15 35rJ‘^°® y'3(4-r2) 

'3r® 2r 9 \ _i —r'^ + r + 3 

'^’.^~15'^y~T5'^^y'^°® ^3(4 -r 2 ) 


Use of Eqs. (C1)-(C3) into Eq. (4.22) gives 


where 


n(l) _ n(l) _|_ n(l) 


(1) 1 ^ 9 2 V , 29g4 41(75 79(76 37g7 

= ( 1 + --9 - 


24 140 560 5040 16 800 2100 25 200 


(C3) 

(C4) 

(C5) 


with 


=- 


22 400g3 


2^2 

T05 




(375 323 -b 1 710 828g -b 548 814g2 _ 204 400 ( 76 ) - tt (23 669 -b 76 404g 


-b56 562g2 -b 12 240g6) -b 135 (387 -b 1164g -b 990g2 -b 272g6) qqs 


-1 


(C6) 


(C7) 


= 


11200g3 


V2-2q- g2 
105 


(375 323 -b 1 599 922g -b 413 153g^ - 465 800g6 - 241 510g^ - 108 524g5 


(I — qU 

-40 954g6 - 4760g^ -b 2555q^ + 1050?® -b 105g^°) -b ^ ^ ' (1603 -b 2243g -b 1479g^ -b 855g6 

-b405g'‘ -b 117g5 -b 17g6 + cos"! -8(1 + gf (1373 _ 336^ _ 72q^ + 96g6 -b 24g'‘) 


V3(l - q) 


_,l + q , (3 + g)5 


X cos-" -b 2 —(111 -b 171g - 117g2 -b 135g3 -b 25g'‘ - 3q^ - Sg® -b g^) cos"! ^ 


V3 


73(3 + g) 


■ (C8) 


Taking into account that (C6) is exact for g < go, and using Eqs. (2.7), (2.8), (2.12), and (3.10), one can obtain the 
following exact expression of E 11112 for g < gp: 


£^11112 = (0^ J (1 + 129 + 66g2 + 40g6 + 81g4 + 


108g5 6516g^ 1521g8 3384g9 6858g 


10 


35 


10 


35 


175 


1458g“ 243gi2\ 27rg3/ ( 1 ) (i)N 

“175-+ 


(C9) 
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